Channel estimation for wireless systems without matrix inversion

ABSTRACT

In various embodiments, techniques are provided to determine channel characteristics of various communication systems such as OFDM systems or systems using a plurality of transmit antennas by using various sets of training symbols that produce zero cross-correlation energy. Channel communication can accordingly be simplified as the zero cross-correlation property allows for channel estimation without a matrix inversion.

This application is a continuation of U.S. patent application Ser. No. 11/508,480 filed on Aug. 22, 2006 now U.S. Pat. No. 7,305,051, entitled “OPTIMUM TRAINING SEQUENCES FOR WIRELESS SYSTEMS” (currently allowed), which is a continuation of U.S. patent application Ser. No. 09/862,755 filed on May 21, 2001, entitled “OPTIMUM TRAINING SEQUENCES FOR WIRELESS SYSTEMS”, which has issued as U.S. Pat. No. 7,103,115. Each of the above cited applications is herein incorporated by reference in their entirety.

BACKGROUND OF THE INVENTION

1. Field of Invention

This invention relates to channel estimation in wireless systems.

2. Description of Related Art

As wireless communications systems are deployed around the world, the importance of providing clear and cost-effective communication services increases. Unfortunately, providing clear communications can require mitigating various obstacles such as inter-symbol-interference (ISI). To reduce ISI, a technique known as orthogonal frequency division multiplexing (OFDM) can be used. Orthogonal frequency division multiplexing is a communication paradigm where a single communication channel is divided into many narrow sub-bands, which then can be transmitted in parallel. By transmitting symbols in this fashion, the duration of each symbol can be dramatically increased, which can greatly reduce or completely eliminate ISI problems.

Unfortunately, individual sub-bands within an OFDM transmission are subject to Rayleigh fading, especially when used in mobile communication systems. While the effects of Rayleigh fading can be mitigated by using multiple transmitter and/or receiver antennas, estimating the channel characteristics for all transmitter-receiver antenna pairs can be difficult and computationally intensive. Accordingly, there is a need for apparatus and techniques that provide for better channel estimation.

SUMMARY OF THE INVENTION

In various embodiments, techniques are provided to determine channel characteristics of various communication systems, such as OFDM systems or systems using a plurality of transmit antennas.

In a first embodiment, methods and apparatus for transmitting training symbols based on a set of first training symbols and one or more sets of second training symbols using one or more second communication channels is provided where a cross-correlation between the first set of training symbols and at least one of the sets of second training symbols is essentially zero.

In a second embodiment, methods and apparatus for channel estimation are provided by first receiving a first set of training symbols related to a first antenna and a set of second training symbols related to a second antenna and then estimating various communication channels. When the cross-correlation between the sets of training symbols is essentially zero, channel estimation can be achieved without a matrix inversion, thus simplifying channel estimation.

BRIEF DESCRIPTION OF THE DRAWINGS

The invention is described in detail with regard to the following figures, wherein like numerals reference like elements, and wherein:

FIG. 1 is a block diagram of an exemplary communication system;

FIG. 2 depicts an OFDM signal having multiple sub-bands;

FIG. 3 depicts an exemplary communication signal of an OFDM sub-band;

FIG. 4 is a block diagram of an exemplary OFDM encoder;

FIG. 5 is a block diagram of an exemplary training symbol generator;

FIG. 6 is a block diagram of an equalizer with an exemplary channel estimator; and

FIG. 7 is a flowchart outlining an exemplary technique for generating and communicating with sets of training symbols.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

This non-provisional application incorporates the substance of U.S. patent application Ser. No. 09/861,811 entitled “CHANNEL ESTIMATION FOR WIRELESS SYSTEMS” to Ye LI now U.S. Pat. No. 7,012,966. The above non-provisional application is concurrently filed and commonly assigned and hereby incorporated by reference in its entirety including all references cited therein.

For wireless systems, channel estimation can be difficult and computationally intensive. Generally, channel estimation for wireless systems can be performed by embedding a pattern of known symbols called training symbols within a transmitted signal and comparing the embedded training symbols against an expected pattern of training symbols at a receiver.

However, for communication systems such as orthogonal frequency division multiplexed (OFDM) systems having multiple sub-bands, each sub-band must be estimated, thus increasing computational requirements. When an OFDM communication system uses multiple transmit antennas and/or multiple receive antennas, each sub-band for each communication channel must be estimated, thus further increasing computational requirements. Fortunately, the problem of estimating the channel parameters for different sub-bands in an OFDM channel can be simplified based on the correlated nature of the sub-bands. Examples of such systems can be found in Yi, L., Seshardi, N. and Ariyavisitakul, S., “Channel estimation for OFDM systems with transmitter diversity in mobile wireless channels”, IEEE Journal of Selected Areas in Comm., Vol. 17, pp. 461-471 (March 1999) incorporated herein by reference in its entirety. Unfortunately, regardless of any advantages posed by intra-channel simplification for a single OFDM channel, the problem of estimating different OFDM channels remains.

However, in systems using multiple transmit antennas, by carefully choosing different patterns of training symbols such that each pattern of training symbols for any transmit antenna has zero cross-correlated energy with the training symbol pattern of any other transmit antenna, channel estimation can be greatly simplified. Such judicious choices of training symbols not only can reduce computational complexity, but also can provide the most accurate channel estimates possible.

FIG. 1 is a block diagram of an exemplary transmission system 100. The transmission system 100 includes an encoder 110 having a number of associated OFDM transmitters 120-1, 120-2, . . . 120-N and respective transmit antennas 130-1, 130-2, . . . 130-N, and an equalizer 160 having a number of associated OFDM receivers 150-1, 150-2, . . . 150-M with respective receive antennas 140-1, 140-2, . . . 140-M.

In operation, the encoder 110 can form blocks of symbols that can be provided to the various OFDM transmitters 120-1, 120-2, . . . 120-N. The various OFDM transmitters 120-1, 120-2, . . . 120-N, in turn, can modulate the blocks of symbols into electromagnetic carriers such as radio-frequency waves that can then be transmitted using their respective transmitting antennas 130-1, 130-2, . . . 130-N. The various radio-frequency signals 135 can then be received by the receive antennas 140-1, 140-2, . . . 140-M and fed to their respective OFDM receivers 150-1, 150-2, . . . 150-M. The OFDM receivers 150-1, 150-2, . . . 150-M can then transform the received radio-frequency signals 135 into base-band signals, digitize the base-band signals and provide the digitized base-band signals to the equalizer 160. The equalizer 160, in turn, can extract the symbols from the digitized base-band signals and perform various operations on the symbols.

As shown in FIG. 1, the radio-frequency signals 135 transmitted by each transmit antenna 130-1, 130-2, . . . 130-N can be subsequently received by each of the receiving antennas 140-1, 140-2, . . . 140-M. While FIG. 1 depicts the various communication channels as single direct paths between each transmit/receive antenna pair, it should be appreciated that each radio-frequency signal 135 can propagate from each transmit antenna 130-1, 130-2, . . . 130-N to each receive antenna 140-1, 140-2, . . . 140-M not only through a direct path, but can also propagate from each transmit antenna 130-1, 130-2, . . . 130-N to each receive antenna 140-1, 140-2, . . . 140-M through a variety of indirect paths (not shown).

The various radio-frequency signal paths for a particular transmit/receive antenna pair can produce a complex communication channel, which can be distinctly different from any other communication channel defined by another transmit/receive antenna pair. Generally, the channel characteristics of an individual mobile wireless channel i.e., the impulse response, can be described by Eq. (1):

$\begin{matrix} {{{h\left( {t,\tau} \right)} = {\sum\limits_{k}\;{{\gamma_{k}(t)}{c\left( {\tau - \tau_{k}} \right)}}}},} & (1) \end{matrix}$ where τ_(k) is the delay of the k-th path, γ_(k)(t) is the corresponding complex amplitude for the k-th path, and c(t) is the shaping pulse for the k-th path whose frequency response is usually a square-root raised cosine Nyquist filter. When a communication channel is a mobile wireless channel, the motion of a vehicle can affect the complex amplitudes γ_(k)(t)'s making each complex amplitude γ_(k)(t) a wide-sense stationary (WSS), narrow-band complex Gaussian process that can be independent for different signal paths.

From Eq. (1), the frequency response H(t,f) of a communication channel at time t can be described by Eq. (2):

$\begin{matrix} {{{H\left( {t,f} \right)}\overset{\Delta}{=}{\int_{- \infty}^{+ \infty}{{h\left( {t,\tau} \right)}{\exp\left( {{- {j2\pi}}\; f\;\tau} \right)}{\mathbb{d}\tau}}}}{or}} & (2) \\ {{\overset{\Delta}{=}{{C(f)}{\sum\limits_{k}\;{{\gamma_{k}(t)}{\exp\left( {{- {j2\pi}}\; f\;\tau_{k}} \right)}}}}}{where}} & (3) \\ {{{C(f)}\overset{\Delta}{=}{\int_{- \infty}^{+ \infty}{{c(\tau)}{{\exp\left( {{- {j2\pi}}\; f\;\tau_{k}} \right)}.}}}}\ } & (4) \end{matrix}$

For OFDM systems with proper cyclic extension and timing, the channel frequency response can be expressed by Eq.(5):

$\begin{matrix} {{{{H\left\lbrack {n,k} \right\rbrack}\overset{\Delta}{=}{{H\left( {{nT}_{f},{k\;\Delta\; f}} \right)} = {\sum\limits_{t = 0}^{K_{o} - 1}\;{{h\left\lbrack {n,l} \right\rbrack}W_{K}^{kl}}}}},}\mspace{11mu}} & (5) \end{matrix}$ where h[n,l]

h(nT_(f), kT_(f)/K), W_(K)=exp(−j2π/K), K is the number of sub-bands (tones) in an OFDM block, T_(f) is the block length and Δf is the sub-band (tone) spacing.

Following Eq. (5), the frequency response at the k-th tone of an n-th block of OFDM symbols corresponding to an i-th transmit antenna can be expressed by Eq.(6):

$\begin{matrix} {{H_{i}\left\lbrack {n,k} \right\rbrack} = {\sum\limits_{l = 0}^{K_{o} - 1}\;{{h_{i}\left\lbrack {n,l} \right\rbrack}{W_{K}^{kl}.}}}} & (6) \end{matrix}$

Equation (6) demonstrates that Hi[n, k] can be obtained by estimating or otherwise acquiring h_(i)[n, k]. Accordingly, the received signal r[n, k] at each receive antenna 140-1, 140-2, . . . 140-M can be expressed by Eq. (7):

$\begin{matrix} {{{r\left\lbrack {n,k} \right\rbrack} = {{\sum\limits_{i = 1}^{M}\;{{H_{i}\left\lbrack {n,k} \right\rbrack}{t_{i}\left\lbrack {n,k} \right\rbrack}}} + {w\left\lbrack {n,k} \right\rbrack}}},} & (7) \end{matrix}$ where M is the number of transmit antennas, k denotes a particular OFDM sub-band and k=0, 1, . . . , K−1 for all n blocks. If the transmitted signals t_(i)[n, k]'s from each transmit antenna contain known signals such as training symbols, the temporal estimation of the various communication channels h_(i)[n, l]'s can be derived using Eq. (8):

$\begin{matrix} {{{\begin{pmatrix} {Q_{11}\lbrack n\rbrack} & {Q_{12}\lbrack n\rbrack} & \cdots & {Q_{1P}\lbrack n\rbrack} \\ {Q_{21}\lbrack n\rbrack} & {Q_{22}\lbrack n\rbrack} & \; & \vdots \\ \vdots & \; & ⋰ & \mspace{11mu} \\ {Q_{P\; 1}\lbrack n\rbrack} & \cdots & \; & {Q_{PP}\lbrack n\rbrack} \end{pmatrix}\begin{pmatrix} {{\overset{\sim}{h}}_{1}\lbrack n\rbrack} \\ {{\overset{\sim}{h}}_{2}\lbrack n\rbrack} \\ \vdots \\ {{\overset{\sim}{h}}_{P}\lbrack n\rbrack} \end{pmatrix}} = \begin{pmatrix} {p_{1}\lbrack n\rbrack} \\ {p_{2}\lbrack n\rbrack} \\ \vdots \\ {p_{P}\lbrack n\rbrack} \end{pmatrix}},{or}} & (8) \\ {{\begin{pmatrix} {{\overset{\sim}{h}}_{1}\lbrack n\rbrack} \\ {{\overset{\sim}{h}}_{2}\lbrack n\rbrack} \\ \vdots \\ {{\overset{\sim}{h}}_{P}\lbrack n\rbrack} \end{pmatrix} = {\begin{pmatrix} {Q_{11}\lbrack n\rbrack} & {Q_{12}\lbrack n\rbrack} & \cdots & {Q_{1P}\lbrack n\rbrack} \\ {Q_{21}\lbrack n\rbrack} & {Q_{22}\lbrack n\rbrack} & \; & \vdots \\ \vdots & \; & ⋰ & \mspace{11mu} \\ {Q_{P\; 1}\lbrack n\rbrack} & \cdots & \; & {Q_{PP}\lbrack n\rbrack} \end{pmatrix}^{- 1}\begin{pmatrix} {p_{1}\lbrack n\rbrack} \\ {p_{2}\lbrack n\rbrack} \\ \vdots \\ {p_{P}\lbrack n\rbrack} \end{pmatrix}}},{where}} & (9) \\ {{{{\overset{\sim}{h}}_{i}\lbrack n\rbrack}\overset{\Delta}{=}\left( {{{{\overset{\sim}{h}}_{i}\left\lbrack {n,0} \right\rbrack}\ldots}\mspace{11mu},{{\overset{\sim}{h}}_{i}\left\lbrack {n,{K_{o} - 1}} \right\rbrack}} \right)^{T}},} & (10) \\ {{{Q_{ij}\lbrack n\rbrack}\overset{\Delta}{=}\left( {q_{ij}\left\lbrack {n,{i - j}} \right\rbrack} \right)_{i,{j = 1}}^{K_{o}}},{and}} & (11) \\ {{{p_{i}\lbrack n\rbrack}\overset{\Delta}{=}\left( {{p_{i}\left\lbrack {n,0} \right\rbrack},\ldots\mspace{11mu},{p_{i}\left\lbrack {n,{K_{o} - 1}} \right\rbrack}} \right)^{T}},} & (12) \end{matrix}$ where {tilde over (h)}_(i)[n] is the estimated channel impulse response for the channel between an i-th transmit antenna and a particular receive antenna, Q_(ij)[n] is a measure of the correlated energy between the set of training symbols t_(i)[n, k] of an i-th transmit antenna and the set of training symbols t_(j)[n, k] of a j-th transmit antenna and p_(i)[n] is the cross-correlation vector between the set of training symbols t_(i)[n, k] of an i-th transmit antenna and a received signal r[n, k] at a particular receive antenna, and where:

$\begin{matrix} {{{q_{ij}\left\lbrack {n,l} \right\rbrack}\overset{\Delta}{=}{\sum\limits_{k = 0}^{K - 1}\;{{t_{i}\left\lbrack {n,k} \right\rbrack}{t_{j}^{*}\left\lbrack {n,k} \right\rbrack}W_{K}^{- {kl}}}}},{and}} & (13) \\ {{p_{i}\left\lbrack {n,l} \right\rbrack}\overset{\Delta}{=}{\sum\limits_{k = 0}^{K - 1}\;{{r\left\lbrack {n,k} \right\rbrack}{t_{i}^{*}\left\lbrack {n,k} \right\rbrack}{W_{K}^{- {kl}}.}}}} & (14) \end{matrix}$

Equation (13) provide a measure of the cross-correlation of energy between two sets of training symbols, i.e., a cross-correlation estimate, and Eq. (14) provides a similar measure between a particular set of training symbols and a received signal. While the particular cross-correlation estimates can be derived using Eqs. (13) and (14), it should be appreciated that any substantially accurate measure of correlation energy derived by any known or later developed means can be used without departing from the spirit and scope of the present invention.

From Eqs. (8) and (9), one can see that a matrix inversion is required to get the temporal estimation of each channel impulse response h_(i)[n, k], i.e., {tilde over (H)}=Q⁻¹ P.

FIG. 2 depicts an exemplary OFDM signal 200 displayed along a time axis 210 and against a frequency axis 220. As shown in FIG. 2, the OFDM signal 200 contains a number of individual sub-bands (tones) 230-1, 230-2, . . . 230-K with each respective sub-band centered on a respective one of closely spaced frequencies f₁, f₂, . . . f_(K). FIG. 3 depicts an exemplary communication signal 300 capable of being embedded in the various sub-bands of FIG. 2. As shown in FIG. 3, the communication signal 300 contains a number of sync symbols 310, a number of training symbols 320, a number of data symbols 330 and a number of guard symbols 340.

Data symbols, also known as payload symbols, can contain information to be transmitted. Guard symbols are symbols that can pad either or both of the beginning and end of a burst transmission and can be used for a variety of purposes including providing buffering, timing and synchronization. Sync symbols are predetermined symbols placed at various strategic positions within a block of data that can allow a receiver to synchronize or otherwise extract timing information from a transmitted signal.

Training symbols, like sync symbols, can be predetermined symbols placed at known positions. However, unlike sync symbols, training symbols are usually configured to enable an equalizer to estimate a given communication channel. It should be appreciated that, in various exemplary embodiments, the training symbols 320 can be any set of symbols suitable for training an equalizer. For example, the exemplary training symbols 320 can be formed taking into account various factors such as their suitability for clock recovery, frequency-shift estimation, their peak-to-average ratio of signal strength or any other known or later recognized factor useful for generating an advantageous or otherwise adequate training sequence.

While the exemplary communication signal 300 is a burst signal with a particular form, it should be appreciated that the form of a burst signal can vary without departing from the spirit and scope of the present invention. For example, it should be appreciated that the training symbols 320 can be dispersed intermittently within the payload symbols 330. It should further be appreciated that while the exemplary communication signal 300 is a burst signal, the communication signal 300 can take various other forms such as a continuous signal in which various training symbols can be periodically embedded.

The exemplary communication signal 300 can be modulated according to a constant modulus scheme such as a quadrature amplitude modulated (QAM) scheme. However, it should be appreciated that the modulation scheme of the communication signal 300 can vary and can take the form of any known or later developed modulation scheme, such as BPSK, QPSK, OPSK, FSK, and the like, without departing from the spirit and scope of the present invention.

FIG. 4 is a block diagram of the exemplary OFDM encoder 110 of FIG. 1 with an associated data source 410. The OFDM encoder 110 includes a sync pattern generator 420, a training symbol generator 430, a guard symbol generator 440 and a number of combining circuits 450.

In operation, the combining circuits 450 can receive a stream of data (payload) symbols from the data source 410 via link 412 and arrange the data stream into a block of K separate sub-streams of data symbols such that each sub-stream of data symbols can be transmitted in a respective OFDM sub-band.

Next, the combining circuits 450 can receive a pattern of guard symbols from the guard symbol generator 440 and append the guard symbols to the end of each data sub-stream to produce a block of K sub-streams of data/guard symbols.

The combining circuits 450 then can replicate the block of data/guard symbol sub-streams N number of times, and retrieve N separate patterns of training symbols {T₁, T₂, . . . T_(N)} from the training symbol generator 430 via link 432. For each of the N blocks, the combining circuits 450 then can replicate each respective pattern of training symbols K number of times and insert the duplicate patterns of training symbols into each of the K separate sub-streams to form N blocks of training/data/guard symbol sub-streams with each of the N blocks having a different pattern of training symbols in each K sub-band.

While the exemplary combining circuits 450 inserts a single pattern of training symbols into each of the K sub-streams for a given block, it should be appreciated that in various embodiments, the training symbol generator 430 can provide K×M different patterns of training symbols such that the combining circuits 450 can assign each sub-stream of training/data/guard symbols a different pattern of training symbols. In still other embodiments, it should be appreciated that each of the N blocks can have a number of patterns of training symbols such that various sub-bands can share one of several available patterns of training symbols. For example, a given OFDM block having thirty-two sub-streams can receive two patterns of training symbols, replicate each pattern of training symbols sixteen times and assign each training symbol pattern to alternating sub-streams.

After each block of training/data/guard symbol sub-streams is formed, the combining circuits 450 can then receive a pattern of sync symbols from the sync pattern generator 420, replicate the sync symbols and append the various training/data/guard symbols to the sync symbols to produce N blocks of K sub-streams of sync/training/data/guard symbols.

While the exemplary OFDM encoder 110 can arrange various symbols to form finite blocks of symbols capable of being transmitted as bursts, it should be appreciated that, in various exemplary embodiments, the OFDM encoder can alternatively form K number of continuous streams of symbols having various training and data symbols, without departing from the spirit and scope of the present invention.

FIG. 5 is a block diagram of the exemplary training symbol generator 430 of FIG. 4. The training symbol generator 430 has a controller 510, a system memory 520, a first training sequence generator 530, a second training sequence generator 540 and an input/output interface 590. The controller 510 interfaces with the various other components 520-590 using control/data bus 502. While the exemplary training symbol generator 430 is depicted as a bussed architecture, it should be appreciated that the functions of the various components 510-590 can be implemented using various other architectures such as complex circuits based on application specific integrated circuits, programmable logic devices, discrete logic and the like.

In operation, and under control of the controller 510, the input/output interface 590 can receive a command to provide N sets of training symbols relating to N separate transmit antennas via link 232 and store the command in the memory 520. In various exemplary embodiments, the input/output interface 590 can receive commands or other information from any device such as a wireless transmitter, a wire transmitter, an optical transmitter, a disc drive, a UART LAN, WAN, parallel digital interface, serial digital interface, software interface or any known or later developed combination of software and hardware without departing from the spirit and scope of the present invention.

After the controller 510 imports and stores the command, the controller 510 can direct the first training sequence generator 530 to generate a first set of training symbols t₁[n, k] according to Eq. (15): t₁[n,k]={S₀,S₁,S₂, . . . S_(K−1)},  (15) where each S_(i) represents a valid symbol state. As discussed above, a set of training symbols such as t₁[n, k] can be any sequence of symbols that is good or otherwise suitable for taking into account factors such as clock recovery, frequency-shift estimation, peak-to-average signal strength ratio or any other known or later recognized factor useful for generating an advantageous or otherwise adequate training sequence.

The first training sequence generator 530 can generate training symbols for communication systems having a constant-modulus signal such that the absolute value of any training sequence is equal to one (|t₁[n, k]|=1). However, it should be appreciated that, in other exemplary embodiments, the first training sequence generator 530 can generate training sequences applicable to any modulation format without departing from the spirit and scope of the present invention.

After the first training sequence generator 530 generates the first set of training symbols t₁[n, k], the controller 510 can provide the first set of training symbols to the second training sequence generator 540. The second training sequence generator 540 can receive the first set of training symbols and generate subsequent sets of training symbols t_(i)[n, k] where i=2, 3, . . . N.

As discussed above, it can be desirable to generate OFDM training sequences such that the cross-correlated energy between any training sequences transmitted by different antennas is essentially zero, i.e., Q_(ij)[n]=0. Accordingly, the second training sequence generator 540 can generate a second training sequence t₂[n, k] such that the cross-correlation energy between t₁[n, k] and t₂[n, k] is essentially zero, i.e., Q₁₂[n]=0. In various exemplary embodiments, such a training sequence t₂[n, k] can be generated by replicating and phase-shifting the first training sequence t₁[n, k] according to Eq. (16):

$\begin{matrix} {{{t_{2}\left\lbrack {n,k} \right\rbrack} = {{{t_{1}\left\lbrack {n,k} \right\rbrack}{\exp\left( {{- j}\frac{2\pi\;{kl}_{0}}{K}} \right)}} = {{t_{1}\left\lbrack {n,k} \right\rbrack}W_{K}^{- {kl}_{o}}}}},} & (16) \end{matrix}$ for some l_(o) with K_(o)≦K−K_(o). Then it can be directly checked that q ₁₂ [n,l]=kδ[l−l _(o)],  (17) where δ[l] denotes the unit impulse function for a frequency l and l₀ is a reference frequency. Eq. (17) implies that q₁₂[u,l]=0 for |l|≦K_(o)−1, and therefore Q₁₂[n]=0 and Q₂₁[n]=Q₁₂ ^(H)[n]=0.

After the second training sequence generator 540 generates the second training sequence t₂[n, k], it can further generate various other training sequences t₃[n, k] t₄[n, k], . . . , t_(N)[n, k] based on the same principle of Eq.(16) such that Q_(ij)[n]=Q_(ji)[n]=0 for all i≠j. Furthermore, like the first training sequence generator 530, the second training sequence generator 540 can generate t_(i)[n, k]'s such that (|t_(i)[n, k]|=1). In general, for communication systems having N(N≦K/K_(o)) transmit antennas, the second training sequence generator 540 can generate the remaining N-1 training sequences according to Eq.(18): t _(i) [n,k]=t ₁ [n,k]W _(K) ^(− K) ^(o) ^((i-1)k),   (18) where

${\overset{\_}{K}}_{o} = {\left\lfloor \frac{K}{N} \right\rfloor \geq K_{o}}$ and └x┘ denotes the largest integer no larger than x. Then for two channels i and j where i<j, the correlation component q_(ij)[n,l] can be described by Eq.(19):

$\begin{matrix} {{{q_{ij}\left\lbrack {n,l} \right\rbrack} = {\sum\limits_{k = 0}^{K - 1}\;{{t_{i}\left\lbrack {n,k} \right\rbrack}{t_{j}^{*}\left\lbrack {n,k} \right\rbrack}W_{K}^{- {kl}}}}},} & (19) \end{matrix}$ and by carefully selecting the relative phases between the training sequences for different transmit antennas, Eq. (9) above can be reduced to the form of Eq. (20) below:

$\begin{matrix} {\begin{pmatrix} {{\overset{\sim}{h}}_{1}\lbrack n\rbrack} \\ {{\overset{\sim}{h}}_{2}\lbrack n\rbrack} \\ \vdots \\ {{\overset{\sim}{h}}_{i}\lbrack n\rbrack} \end{pmatrix} = {\begin{pmatrix} {Q_{11}\lbrack 0\rbrack} & 0 & \ldots & 0 \\ 0 & {Q_{22}\lbrack 1\rbrack} & \; & 0 \\ \vdots & \; & ⋰ & \vdots \\ 0 & 0 & \ldots & {Q_{ij}\lbrack n\rbrack} \end{pmatrix}^{- 1}\begin{pmatrix} {p_{1}\lbrack n\rbrack} \\ {p_{2}\lbrack n\rbrack} \\ \vdots \\ {p_{i}\lbrack n\rbrack} \end{pmatrix}}} & (20) \end{matrix}$

Because all off-axis elements of the Q matrix are essentially zero, determining the inverse matrix Q⁻¹ can be greatly simplified. Furthermore, for communication systems using a constant modulus modulation, i.e., |t_(i)[n, k]|=1, then Q in Eq. (20) can be reduced to Q=K×I, for all i=1, 2, . . . N, and Eq. (20) can be further reduced to the form of Eq. (21):

$\begin{matrix} {{\begin{pmatrix} {{\overset{\sim}{h}}_{1}\lbrack n\rbrack} \\ {{\overset{\sim}{h}}_{2}\lbrack n\rbrack} \\ \vdots \\ {{\overset{\sim}{h}}_{i}\lbrack n\rbrack} \end{pmatrix} = {\begin{pmatrix} K & 0 & \ldots & 0 \\ 0 & K & \; & \vdots \\ \vdots & \; & ⋰ & 0 \\ 0 & \ldots & 0 & K \end{pmatrix}^{- 1}\begin{pmatrix} {p_{1}\lbrack n\rbrack} \\ {p_{2}\lbrack n\rbrack} \\ \vdots \\ {p_{i}\lbrack n\rbrack} \end{pmatrix}}},\mspace{14mu}{or}} & (21) \\ {\begin{pmatrix} {{\overset{\sim}{h}}_{1}\lbrack n\rbrack} \\ {{\overset{\sim}{h}}_{2}\lbrack n\rbrack} \\ \vdots \\ {{\overset{\sim}{h}}_{i}\lbrack n\rbrack} \end{pmatrix} = {\frac{1}{K}\begin{pmatrix} 1 & 0 & \ldots & 0 \\ 0 & 1 & \; & \vdots \\ \vdots & \; & ⋰ & 0 \\ 0 & \ldots & 0 & 1 \end{pmatrix}\begin{pmatrix} {p_{1}\lbrack n\rbrack} \\ {p_{2}\lbrack n\rbrack} \\ \vdots \\ {p_{i}\lbrack n\rbrack} \end{pmatrix}}} & (22) \end{matrix}$

Accordingly, the problem of determining the set of channel characteristics H_(i)[n,l] for a communication system having N transmit antennas can be reduced to solving Eq. (23):

$\begin{matrix} {{{h_{i}\left\lbrack {n,l} \right\rbrack} = {\frac{1}{K}{p_{i}\left\lbrack {n,l} \right\rbrack}}},{{or}\mspace{14mu}{alternatively}}} & (23) \\ {{h_{i}\left\lbrack {n,l} \right\rbrack} = {\frac{1}{K}{p_{i}\left\lbrack {n,{l - {\left( {i - 1} \right){\overset{\_}{K}}_{o}}}} \right\rbrack}}} & (24) \end{matrix}$ for l=0, . . . , K _(o)−1 and i=1, 2 . . . , N. Therefore, by carefully selecting the relative phases between the training sequences for different transmit antennas, the timing sequences for each channel of each respective transmit antenna on p_(i)[n, l] are shifted to different regions in the time domain and the parameters for different channels can be easily estimated without using a computationally intensive matrix inversion.

As discussed above, for systems where every off-axis element of Q is exactly zero, i.e., Q_(ij)[n]=Q_(ji)[n]=0 for all i≠j, it follows that Q⁻¹ can be directly replaced by (1/K) I. However, if one or more off-axis Q_(ij)[n] are close to, but not exactly, zero, Q⁻¹ may still be approximated by (1/K) I. That is, it should be appreciated that, in various embodiments, Q_(ij)[n] can be any value of Q_(ij)[n]=Q_(ji)[n]≈0 such that Q⁻¹ containing such Q_(ij)[n] can still be functionally replaced with (1/K) I, with an understanding that some performance degradation may occur.

A further advantage to designing the various sets of training symbols according to the above-described technique is that, when Q_(ij)[n]=0 for i≠j as in Eqs. (21) or (22), a channel estimator estimating the various channels can effectively attain the theoretical lower mean square error (MSE) boundary. Accordingly, such sets of training symbols can also enable an estimator to achieve the best theoretical performance.

While the above-described training symbol design technique can be applied to various OFDM systems, it should be appreciated that the above-identified technique can be applied to any other communication system using multiple transmit antennas. Furthermore, the training sequences developed above can be easily adapted to pilot sequence design for pilot symbol aided channel estimation or channel estimation in single carrier systems. For pilot symbol aided channel estimation in OFDM systems with pilot tones scattered into different times and frequencies, the pilot sequences can be described as two-dimensional and the optimum sequence design strategy can accordingly be directly used. Furthermore, the above-described approach can be more flexible than conventional approaches to pilot symbol aided channel estimation since the relative phases of the pilot sequences for different transmit antennas can be shifted in a two-dimensional claim.

FIG. 6 is a block diagram of the exemplary equalizer 160 of FIG. 1. The exemplary equalizer 160 includes a diversity gain processor 610 and a channel estimator 620 containing a training symbol database 630.

In operation, the diversity gain processor 610 and channel estimator 620 can receive various received signals r₁[n, k], r₂[n, k], . . . r_(M)[n, k] such as blocks of symbols from OFDM transmissions via links 152-1, 152-2, . . . 152-M. As discussed above, each of the received signals r_(i)[n, k] can contain multiple transmit signals t_(i)[n, k] according to Eq. (7) above transmitted by a plurality of transmit antennas, with each transmit signal t_(i)[n, k] having an assortment of sync symbols, training symbols, data symbols and guard symbols, as well as any other symbol types. Also as discussed above, the training symbols embedded in the communication signals for each of the transmit signals t_(i)[n, k] can be known symbol patterns and formed according to Eqs. (15)-(18) above such that there is essentially no cross-correlated energy between each set of training symbols, i.e. Q_(ij)[n, l]=0 for all i≠j.

During operation, channel estimator 620 can extract various sets of expected patterns of training symbols t_(i)[n, k] from the training symbol database 630 and process each received signal r_(i)[n, k] based on the technique described above in Eq. (14) above to derive p_(i)[n, l], which can then be used to determine the various channel characteristics based on Eq. (18) or, in the case where a constant modulus modulation scheme is used, based on Eqs. (21)-(22).

After the channel estimator 620 determines the appropriate set of channel characteristics, the channel estimator 620 can export the set of channel characteristics to the diversity gain processor 610, which can use the set of channel characteristics to provide spatial and temporal equalization for the received signals to produce an estimated symbol stream that can be exported to an external device (not shown) via link 162.

The exemplary diversity gain processor 610 uses a minimum-mean-square-error (MMSE) technique to provide equalization. However, it should be appreciated that the diversity gain processor 610 can use any known or later developed technique useful for temporal and/or spatial equalization without departing from the spirit and scope of the present invention.

FIG. 7 is a flowchart outlining a first exemplary operation for developing and using optimum training sequences. The process starts in step 710 where a first set of training symbols is developed. As described above, the first set of training symbols can be any set of training symbols that is advantageous or otherwise good for characterizing a communication channel. Next, in step 720, a determination is made as to whether to generate another set of training symbols. If another set of training symbols is to be generated, control jumps to step 770; otherwise, control continues to step 730.

In step 770, a next set of training symbols is generated based on the first set of training symbols such that the cross-correlation energy between the first and second set of training symbols is essentially zero. Control then jumps back to step 720 where another determination is made as to whether to generate another set of training symbols. If more sets of training symbols are to be generated, step 770 is repeated such that each new set of training symbols generated will not appreciably have any cross-correlation energy with any of the previously generated sets of training symbols. The exemplary technique generates subsequent timing sequences by phase shifting the first time sequence according to Eqs. (13)-(16) above. However, it should be appreciated that, in various exemplary embodiments, the particular technique used to generate subsequent sets of training symbols can vary without departing from the spirit and scope of the present invention.

In step 730, because no further sets of training symbols are to be generated, the training sequences generated in step 710, 720 and 770 are transmitted. The exemplary sets of training symbols are transmitted according to an OFDM paradigm with the different sets of training symbols transmitted using different transmit antennas. However, as discussed above, it should be appreciated that the various sets of training symbols can be transmitted according to any other known or later developed paradigm using multiple transmit antennas and/or sub-bands without departing from the spirit and scope of the present invention.

Next, in step 740, the training symbols are received. While the exemplary technique uses multiple OFDM receivers coupled to an equalizer, it should be appreciated that the number of OFDM receivers can vary without departing from the spirit and scope of the present invention. Furthermore, as discussed above, while the exemplary receiver can operate according to an OFDM paradigm, it should be appreciated that the receiver can operate according to any other known or later developed paradigm without departing from the spirit and scope of the present invention.

Then, in step 750 the communication channels between the various transmit antennas and receive antennas can be estimated. As discussed above, the exemplary technique can estimate the various channels using a correlation matrix Q having all of the off-diagonal elements equal to zero to make the process of channel estimation extremely simple. Furthermore, as discussed above, in systems using a constant modulus modulation approach, the correlation matrix Q can be reduced to the identity matrix, thereby alleviating the necessity for performing an inverse transform altogether. Control then continues to step 760 where the operation stops.

As shown in FIG. 1-6, the systems and methods of this invention are preferably implemented on a digital signal processor (DSP) or other integrated circuits. However, the systems and methods can also be implemented using any combination of one or more general purpose computers, special purpose computers, program microprocessors or microcontroller and peripheral integrating circuit elements, hardware electronic or logic circuits such as application specific integrated circuits (ASICs), discrete element circuits, programmable logic devices such as PODs, POAs, FPGAs, PALs, or the like. In general, any device on which exists a finite state machine capable of implementing the various elements of FIGS. 1-6 and the flowchart of FIG. 7 can be used to implement the training sequence functions.

While this invention has been described in conjunction with the specific embodiments thereof, it is evident that many alternatives, modifications and variations will be apparent to those skilled in the art. Accordingly, preferred embodiments of the invention as set forth herein are intended to be illustrative, not limiting. There are changes that may be made without departing from the spirit and scope of the present invention. 

1. A method for communicating, comprising: transmitting via a first transmitter a set of first training symbols using a first communication channel; and transmitting via a second transmitter one or more sets of second training symbols using one or more second communication channels; wherein the one or more sets of second training symbols are based on the set of first training symbols and a cross-correlation estimate between the first set of training symbols and at least one of the one or more sets of second training symbols is essentially zero, where a channel estimation is performed without performing a matrix inversion.
 2. The method of claim 1, wherein at least the first set of training symbols is transmitted using an orthogonal frequency division multiplexing technique.
 3. The method of claim 1, wherein each cross-correlation estimate between the first set of training symbols and every set of the one or more sets of second training symbols is essentially zero.
 4. The method of claim 3, wherein each cross-correlation estimate between every two sets of training symbols of the one or more sets of second training symbols is essentially zero.
 5. The method of claim 1, wherein at least one set of the one or more sets of second training symbols is substantially identical to the set of first training symbols with a phase shift.
 6. The method of claim 5, wherein every set of the one or more sets of second training symbols is substantially identical to the set of first training symbols with a respective phase shift.
 7. The method of claim 1, wherein the first set of training symbols is transmitted using a first transmit device and at least one of the one or more sets of second training symbols is transmitted using a second transmitting device.
 8. The method of claim 7, further comprising: receiving the set of first training symbols; receiving at least one of the one or more sets of second training symbols; and characterizing two or more communication channels based on the set of first training symbols and the one or more second sets of training symbols.
 9. The method of claim 8, wherein said characterizing the two or more communication channels does not use a matrix inversion.
 10. A method of communicating, comprising: receiving via a first receiver a set of first training symbols; receiving via a second receiver one or more sets of second training symbols; and characterizing two or more communication channels based on the set of first training symbols and the one or more sets of second training symbols; wherein a cross-correlation estimate between the set of first training symbols and at least one of the sets of second training symbols is essentially zero, where a channel estimation is performed without performing a matrix inversion.
 11. The method of claim 10, wherein at least the set of first training symbols transmitted using an orthogonal frequency division multiplexing technique.
 12. The method of claim 11, wherein each cross-correlation estimate between the set of first training symbols and every set of the one or more sets of second training symbols is essentially zero.
 13. The method of claim 12, wherein each cross-correlation estimate between every two sets of training symbols of the one or more sets of second training symbols is essentially zero.
 14. The method of claim 13, wherein at least one set of the one or more sets of second training symbols is substantially identical to the set of first training symbols with a phase shift.
 15. The method of claim 14, wherein every set of the one or more sets of second training symbols is substantially identical to the set of first training symbols with a respective phase shift.
 16. The method of claim 10, wherein said characterizing two or more communication channels does not use a matrix inversion.
 17. The method of claim 10, wherein the set of first training symbols is transmitted using a first transmit device and the one of the one or more sets of second training symbols is transmitted using a second transmitting device.
 18. An apparatus for communicating, comprising: a first transmit device having a first transmit antenna that transmits a set of first training symbols; and a second transmit device having a second transmit antenna that transmits a set of second training symbols; wherein a cross-correlation estimate between the set of first training symbols and the set of second training symbols is essentially zero, where a channel estimation is performed without performing a matrix inversion.
 19. The apparatus of claim 18, wherein the set of second training symbols is substantially identical to the set of first training symbols with a phase shift.
 20. An apparatus for communicating, comprising: a receive device having at least one receive antenna that receives at least a set of first training symbols transmitted by a first transmit device having a first transmit antenna and a set of second training symbols transmitted by a second transmit device having a second transmit antenna; and an estimator that estimates at least a first channel related to the first transmit device based on at least the set of first training symbols; wherein a cross-correlation estimate between the set of first training symbols and the set of second training symbols is essentially zero, where a channel estimation is performed without performing a matrix inversion.
 21. The apparatus of claim 20, wherein the estimator further estimates the first channel based on at least the set of second training symbols.
 22. The apparatus of claim 21, wherein the estimator estimates the first channel without using said matrix inversion.
 23. The apparatus of claim 22, wherein the set of second training symbols is substantially identical to the set of first training symbols with a phase shift. 